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Spin-dependent structure functions in nuclear matter and the polarized EMC effect 
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An excellent description of both spin-independent and spin-dependent quark distributions and 
structure functions has been obtained with a modified Nambu— Jona-Lasinio model, which is free 
of unphysical thresholds for nucleon decay into quarks — hence incorporating an important aspect 
of confinement. We utilize this model to investigate nuclear medium modifications to structure 
functions and find that we are readily able to reproduce both nuclear matter saturation and the 
experimental F^ N / F2N ratio, that is, the EMC effect. Applying this framework to determine gf p , 
we find that the ratio gt v / g\ v differs significantly from 1, with the quenching caused by the nuclear 
medium being about twice that of the spin-independent case. This represents an exciting result, 
which if confirmed experimentally, will reveal much about the quark structure of nuclear matter. 

PACS numbers: 25.30Mr, 13.60Hb, 24.85+p, 11.80Jy, 12.39Fe, 12.39Ki 



The discovery in the early 80's by the European Muon 
Collaboration (EMC) that nuclear structure functions 
differ substantially from those of free nucleons [1-3] 
caused a shock in the nuclear community. Despite many 
attempts to understand this effect in terms of binding 
corrections it has become clear that one cannot under- 
stand it without a change in the structure of the nucleon- 
like quark clusters in matter [4-6]. Mean- held models of 
nuclear structure built at the quark level, which have 
been developed over the past 15 years, are yielding a 
quantitative description of the EMC effect. Most recently 
it has been demonstrated that at least one of these mod- 
els leads naturally to a Skyrme-type force, with param- 
eters in agreement with those found phenomenologically 
to describe a vast amount of nuclear data [7] . 

A second major discovery by the EMC concerned the 
so-called "spin crisis" [8], which corresponds to the dis- 
covery that the fraction of the spin of the proton car- 
ried by its quarks is unexpectedly small. This has led to 
major new insights into the famous U(l) axial anomaly, 
prompting many new experiments. With this back- 
ground, it is astonishing that, in the 17 years since the 
discovery of the spin crisis, there has been no experimen- 
tal investigation of the spin-dependent structure func- 
tions of atomic nuclei. Of course, such experiments are 
more difficult because the nuclear spin is usually carried 
by just a single nucleon and hence the spin dependence 
is an 0(1/ A) effect. Nevertheless, as we shall see, such 
measurements promise another major surprise, with at 
least one model — which reproduces the EMC effect in 
nuclear matter — suggesting a modification of the spin 
structure function of a bound proton in nuclear matter 
roughly twice as large as the change the spin-independent 
structure function. 

Models of nuclear structure like the quark meson cou- 
pling (QMC) model, achieve saturation through the self- 
consistent change in the quark structure of the color- 



less, nucleon-like constituents — in particular, through 
its scalar polarizability [7, 9]. Physically the idea is ex- 
tremely simple, light quarks respond rapidly to oppose 
an applied scalar field. Specifically, the lower compo- 
nents of the valence quark wave functions are enhanced 
and this in turn reduces the effective aN coupling. The 
fact that changes in the structure of bound nucleons are 
so difficult to find appears to be a result of this mecha- 
nism being extremely efficient and hence yielding only a 
small change in the dominant upper components of the 
valence quark wave functions. 

On the other hand, the spin structure functions are 
particularly sensitive to the lower components and this 
is why the measurement of the spin-dependent EMC ef- 
fect is so promising. Our calculations are made within 
the framework developed by Bentz, Thomas and collab- 
orators [10, 11], in which proper-time regularization [12- 
14] is applied to the NJL model in order to simulate the 
effects of confinement. This model exhibits similar prop- 
erties to the QMC model with the advantage that it is 
covariant. Once we include both scalar and axial- vector 
diquarks, it readily describes nuclear saturation at the 
correct energy and density. Moreover it yields PDFs for 
the free nucleon [15] which are in excellent agreement 
with existing experimental data. 

We write the spin-dependent light-cone quark distri- 
bution of a nucleus with mass number A and helicity H 
as the convolution 



(*a) = 
dy A dx 6(x A 



-y A x)Af q/N (x) Af ( N H / A (y A ), (1) 



where Af q / N (x) is the spin-dependent quark light- 
cone momentum distribution in the bound nucleon, 
Afpf, A (y A ) the light-cone momentum distribution of the 
nucleon in the nucleus and x A € [0, A] is the Bjorken scal- 
ing variable for the nucleus. There have been numerous 



investigations of Af^A (y^) [16] and it is straightfor- 
ward to calculate for any particular nucleus. Examples 
of greatest experimental interest would be single proton 
particle or hole states like 7 Li, n B and 15 JV. In this 
analysis, as our primary focus is the change in Af q /^ in- 
medium, we incorporate the Fermi motion effects on the 
bound proton structure function by replacing Af N ^ A (i/A) 
with the spin-independent distribution Jm/a^a), calcu- 
lated in infinite nuclear matter [11]. 

To calculate Af q / N (x) in our model, it is convenient 
to express it in the form [17, 18] 

Af g/N (x) = -iJ 

Tr( 7 + 75 M(p,fc)), (2) 

where M{p, k) is the quark two-point function in the nu- 
cleoli. Within any model that describes the nucleon as 
a bound state of quarks, this distribution function can 
be associated with a straightforward Feynman diagram 
calculation, where the propagators include the self con- 
sistent scalar and vector fields in the nucleus. 

It is demonstrated in Ref. [11] that the in-medium 
changes to a free nucleon quark distribution can be in- 
cluded as follows. The effect of the scalar field is in- 
corporated by simply replacing the free masses with the 
effective masses in nuclear medium, giving the distribu- 
tion Af q /N (x) [22] and the Fermi motion of the nucleon 
is included by convoluting this distribution (Af q /j^ (x)) 
with the Fermi smearing function, fN/Ao{yA), producing 
the distribution 

^f q /Ao{xA) = J dy A J dz 

5 (x A - y A z) Af q/NQ (z) f N /Ao{VA)- (3) 

The effect of the vector field is then incorporated via the 
scale transformation 

Af q /A (x A ) = |r; M q /A0 (xa = ~jjr x A - , (4) 

where £p = \]yp + Mjtf + 3T / o = Ep + 3Vb is the Fermi 
energy of the nucleon, pp the Fermi momentum and Vq is 
the zeroth component of the vector field felt by a quark. 

To calculate the spin-dependent quark distribution in 
the nucleon, Af q / N0 (x), we use the NJL model to de- 
scribe the nucleon as a quark-diquark bound state, taking 
into account both scalar (j 71 ^ = + ,T = 0, colour 3) and 
axial- vector ( J 71 " = 1 + , T = 1, colour3 ) diquark chan- 
nels. Details of these free space calculations, along with 
a description of the proper-time regularization scheme 
used throughout this paper, may be found in Ref. [15]. 
In short, the quark distribution functions are determined 
from the Feynman diagrams of Fig. 1, with the resulting 
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FIG. 1: Feynman diagrams representing the spin-dependent 
quark distributions in the nucleon, needed to determine 
Af q / N (x), given in Eq. (2). The single line represents 
the quark propagator and the double line the diquark 
i-matrix. The shaded oval denotes the quark-diquark 
vertex function and the operator insertion has the form 

7 + 75 s(x ~ jr~) \ (1 i T z)- The second diagram, which we 
refer to as the "diquark diagram", symbolically represents 
two diagrams, each with the operator insertion on a different 
quark line within the diquark. 

distribution Af q / N0 (x), having no support for negative 
x. Hence, this is essentially a valence quark picture. 

By calculating these Feynman diagrams using the ef- 
fective (density dependent) masses obtained from the nu- 
clear matter equation of state (discussed below) and per- 
forming the transformation, Eq. (4), to include the mean 
vector field, we obtain the spin-dependent u and d dis- 
tributions in a bound proton. Separating the isospin fac- 
tors, gives 

Au A (x) = Af s q/N (x) + l - Af q{D)/N {x) + i Af a q/N {x) 

+ I A fq(D )/N (x) + ^=Af% D)/rf {x), (5) 

Ad A (x) = \ Ar q(D)/N {x) + 2 - AQ N (x) 

+ I A fq(D )/N (x) ~ ^=Af% D)/ff (x). (6) 

The superscripts s, a and m refer to the scalar, axial- 
vector and mixing terms, respectively, the subscript q/N 
implies a quark diagram and q(D)/N a diquark dia- 
gram. Because the scalar diquark has spin zero, we have 
Af^ D y N (x) — and hence the polarization of the d 
quark arises exclusively from the axial-vector and the 
mixing terms. 

The NJL model is a chiral effective quark theory that 
is characterized by a 4-Fermi contact interaction. Us- 
ing Fierz transformations any 4-Fermi interaction can be 
decomposed into various interacting qq and qq channels 
[19]. The terms relevant to this discussion are 

C = ip (i $ — m) ip 

+ ((V^f - (^75 r^) 2 ) - G w (^ 7 M V-) 2 + • ■ ■ 

+ G s (? 75 Cr 2 /3 A ^ T ) (V^C 1 7572/3^) 

+ G a (^^Ct iT2 I3 a ^ (i> T C- l YT 2 Ti(3 A Tlj) , (7) 
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FIG. 2: Spin-dependent quark distributions, Au v and Ad v , 
at the model scale, Qq — 0.16 GeV 2 . There are four curves for 
each quark flavour, with the positive curves representing the 
up distributions. The dotted line is the free nucleon distribu- 
tion, the dot-dash line illustrates the effect of replacing the 
free masses with the effective ones. This distribution convo- 
luted with the Fermi smearing function, Eq. (3), is presented 
as the dashed line, and the final result where the vector field 
is also included via the scale transformation, Eq. (4), is rep- 
resented by the solid line. 



where m is the current quark mass, = y | A (A = 

2,5,7) are the colour 3 matrices and C = 17270- In the 
qq channel we include scalar, pseudoscalar and vector 
components and in the qq channel we have the scalar 
and axial-vector diquarks. The scalar qq interaction term 
generates the scalar field, that is, the constituent quark 
mass M (vacuum value Mo) via the gap equation. The 
vector qq interaction will be used to generate the vector 
field in-medium. The qq interaction terms give the di- 
quark ^-matrices whose poles correspond to the masses 
of the scalar and axial- vector diquarks. The nucleon ver- 
tex function and mass, Mjf, are obtained by solving the 
homogeneous Faddeev equation for a quark and a diquark 
[15]. Because we need to solve this equation many times 
to obtain self-consistency, we approximate the quark ex- 
change kernel by a momentum independent form (static 
approximation) . This necessitates the introduction of an 
additional parameter, c, as explained in Ref. [10]. 

To calculate the mean scalar and vector fields, we need 
the equation of state for nuclear matter. This can be rig- 
orously derived for any NJL Lagrangian using hadroniza- 
tion techniques, but in a simple mean-field approximation 
the result for the energy density has the following form 
[10]: 



£ = £v - 



V 2 
4G W 



d 3 p 
(2^F 



Q(p F -\p\)e p , (8) 



where e p — \J p 2 + Mfj + 3Vb and the vacuum term £y 
has the familiar "Mexican hat" shape. 

The parameters of the model are Air, Ajjy, -Mo, c, G^, 




FIG. 3: Ratios of the spin-independent and spin-dependent 
nuclear to nucleon structure functions at nuclear matter den- 
sity. The top curve is the usual EMC ratio F 2 a n /F2n, where 
F2N is the isoscalar structure function and the superscript A 
represents the in-medium result. The EMC data for nuclear 
matter is taken from Ref. [21]. Our prediction for the polar- 
ized EMC effect, gt p /gi v , is the lower curve. Clearly we find 
a significant effect. 



G s , G a and G u , where Am and Ajjy are the infrared and 
ultraviolet cutoffs used in the proper-time rcgularization. 
The infrared scale is expected to be of the order Aqcd 
and we set it to A//? = 0.28 GeV. We also choose the free 
constituent quark mass to be Mo = 400 MeV [23] and use 
this constraint to fix the static parameter, c. The remain- 
ing six parameters are fixed by requiring = 93 MeV, 
tti-jy = 140 MeV, M N = 940 MeV, the saturation point 
of nuclear matter (p B ,E B ) = (0.17fm~ 3 , 15.7 MeV) and 
lastly the Bjorken sum rule at zero density to be sat- 
isfied, with g A = 1.267. We obtain A uv = 0.66 GeV, 
c = 0.95 GeV, = 17.81 GeV~ 2 , G s = 8.41 GeV~ 2 , 
G a = 1.36 GeV" 2 and G w = 5.58 GeV~ 2 . 

With these model parameters the diquark masses at 
zero density are M s = 0.65 GeV and M a = 1.2 GeV 
and vector field strength is Vo = 0.044 GeV. At satura- 
tion density the effective masses become M* = 0.32 GeV, 
A/ s * = 0.52 GeV, M* = 1.1 GeV and M* N = 0.75 GeV. 

The results for the u and d spin-dependent quark dis- 
tributions, at the model scale, are presented in Fig. 2. 
There are four curves for each quark flavour, represent- 
ing the different stages leading to the full nuclear matter 
result. 

Using these quark distributions we are able to con- 
struct the structure functions, g\ p and gf p , where 
the superscript A represents a structure function in 
the nuclear medium. Analogous results for the spin- 
independent quark distributions [15] allow us to deter- 
mine the isoscalar structure functions F2N and F^ N , and 
hence determine the EMC effect. Evolving [20] these 
distributions to a scale of 10 GeV 2 , we give in Fig. 3 
our results for the ratios F^ N /F2N and g^ p /gi P , that is 
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FIG. 4: Ratio of the quark distributions in nuclear mat- 
ter to the corresponding free distributions, at a scale of 
Q 2 = lOGeV 2 . The solid line represents Au A (x)/Au(x) and 
the dot-dashed line Ad A (x) / Ad(x) . Note, these distributions 
are the full quark distributions and hence include anti-quarks 
generated through Q 2 evolution. 



the EMC and the polarized EMC effect. In the valence 
quark region, the model is able to reproduce the spin- 
independent EMC data extremely well. For the polarized 
ratio we find a significant effect, of the order twice the 
size of the unpolarized EMC effect. 

The nuclear quenching effects on the individual quark 
flavours is presented in Fig. 4. We find that the effect on 
both the u and d distributions is large and approximately 
equal over the valence quark region. The resemblance 
between g^ p /gi p and the ratio Au A (x)/Au(x) is simply 
because the up distribution is enhanced by a factor 4 
relative to the down and strange distributions in proton 
structure functions. Absent from our model is the U(l) 
axial anomaly and sea quarks (at the model scale), which 
prevents a reliable description of structure functions at 
low x. For this reason in Figs. 3 and 4 we do not plot 
our results in this region. 

A thorough understanding of how nuclear medium ef- 
fects arise from the fundamental degrees of freedom — 
the quarks and gluons — represents an important chal- 
lenge for the nuclear physics community. An experimen- 
tal measurement of the polarized EMC effect would be 
another important step toward this goal, providing im- 
portant insights into the quark polarization degrees of 
freedom within a nucleus. Our prediction of a remark- 
ably large signature suggests that this measurement is 
feasible, and if these results are confirmed experimentally 
would yield vital, new information on quark dynamics in 
the nuclear medium. 
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Our results do not depend strongly on this choice, re- 
maining almost unchanged with Mo is between 350 and 
450 MeV. 



